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Abstract 

We study charged lepton flavor violating processes in the context of supersymmetric mod- 
els with non-universal gaugino masses. We consider the generic case of free gaugino masses, 
as well as specific models based on gaugino mediation and SU{^) GUTs. Under the assump- 
tion that the main source of lepton flavor violation resides in the off-diagonal elements of the 
left handed slepton mass matrix, we show that BR(/x — > 67) strongly depends on gaugino 
non-universality. When Mi and M2 have opposite signs, cancellations occurring between 
the neutralino and the chargino contributions can suppress, even completely, the branching 
ratio of charged lepton flavor violating decays. Moreover, we find that for M3 > M2,Mi, 
BR(/i — 67) is an increasing function of M^. 



1 Introduction 



Universal gaugino masses are not a consequence of the supergravity framework but rather an 
additional assumption. In supergravity models, supersymmetry is broken in a hidden sector 
which consists of fields which couple to the MSSM fields only gravitationally. In the minimal 
scenario (mSUGRA), it is assumed that the vacuum expectation value of the gauge kinetic 
function does not break the unifying gauge symmetry. This simplification leads to universal 
gaugino masses. In mSUGRA it is also assumed that all scalars acquire a common mass, 
and that the trilinear couplings are universal. Such simplifying assumptions, however, are not 
necessarily valid in specific scenarios. Indeed, in the literature there is plenty of models which do 
not fulfill such requirements [J. It is then worth to consider the phenomenological implications 
of non minimal scenarios, and in particular of gaugino non- universality [2]. 

So far, the discussion of models with non-universal gaugino masses has focused on the im- 
plications for SUSY particles detection in accelerator experiments (S] and dark matter [1]. One 
of the features of these models is that the cosmological upper bound on the relic density [Sj 
can be easily fulfilled. In fact, the lightest supersymmetric particle may be a wino- or higgsino- 
dominated rather than a bino-dominated neutralino. In both cases, annihilations into W~^W~ 
and Z^Z^, as well as coannihilations with the lightest chargino and (for the higgsino) with 
the next-to-lightest neutralino can suppress the relic density below the currently estimated dark 
matter content of the Universe [SllZj. Additionally, such neutralinos, due to their large couplings 
to the W and Z bosons, typically yield large direct and indirect detection rates, and therefore 
constitute appealing candidates for WIMP search experiments (HI El El- Lepton flavor violating 
(LFV) decays, on the other hand, have never been studied in this context. 

The process /U — > 67 is an ideal candidate to look for new physics jHllSl- Recent experiments 
on neutrino oscillations have shown that flavor is not conserved in the lepton sector, implying 
that lepton flavor violating decays such as — > 67, r — > //7 and /i — > 3e are in principle allowed. 
Since in standard model-like theories these processes are extremely suppressed, their detection 
would be an unmistakable signal of physics beyond the Standard Model. In a supersymmetric 
framework, the precise rate of LFV decays depends on the mechanism that communicates the 
violation of lepton flavor to the soft breaking Lagrangian. To maintain our discussion as general 
as possible, instead of studying one of the different scenarios considered in the literature, we will 
conveniently parametrize the terms assumed to be responsible for LFV at low energies. 

In this article we study the effects that non universal gaugino masses can have in the decay 
// — > 67. In the next section we explain our assumptions and introduce some notation. Then, 
we briefly review the predictions for BR(;U 67) in models with universal gaugino masses, and 
afterwards we proceed to consider general models with non universality, which we study in some 
detail. Finally, we discuss models based on the unifying group SU{5), which predict a fixed 
pattern of non-universality, and models with gaugino mediated SUSY breaking. 

2 Generalities 

A generic model with non universal gaugino masses contains, in addition to the mSUGRA param- 
eters mo,ao,tan/3 and sign(/x), three gaugino masses Mi, M2, corresponding, respectively, to 
the C/(l), SU{2) and SU{3) gauge groups. The effect of non universality can be parametrized, 
at the high energy GUT scale Mqut ^ 2 10^'^ GeV, by the dimensionless ratios r2 = M2/M1 
and r3 = M3/M1. In the following, we will study the dependence of BR{^ — > 67) on r2 and 
for different values of the remaining parameters. 

At low energies, the only possible sources of lepton fiavor violation in the MSSM soft breaking 
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Lagrangian are non- vanishing off-diagonal elements in m|, m| and Ag. We will assume that 
the dominant source of LF violation resides in the left-handed slepton mass matrix m|, and 
parametrize the off-diagonal element (m|)2i in terms of the diagonal entry as 

("l|)21 = <5od("i|)i1 = '^OD("i|)22, (1) 

600 being a free parameter. The off-diagonal elements (m|)i3 and (rn-|)23 play a subdominant 
role in the process fi — > 67, since they only contribute through double mass insertions, and 
will therefore be disregarded in what follows. The matrices m| and A^,, on the other hand, are 
taken to be diagonal. The outlined scenario applies, for instance, in see-saw inspired models 
of LFV im in]. We remark that, though our results refer to the decay fi — > ej, with obvious 
substitutions, they also hold for the decays r — > /i7 and r — > 67. 
The amplitude for the process /x — > 67 is written 

u^ip) (2) 

where q is the momentum of the photon, e is the electric charge, e* the photon polarization 
vector. Me and are the wave functions for the electron and the muon, and p is the momentum 
of the muon. Each coefficient in the above equation can be written as a sum of two terms 

where ^(")^'^ and A^'^^^'^ stand, respectively, for the contribution from the neutralino-slepton 
and the chargino-sneutrino loops. Explicit expressions for these contributions can be found in 
jlUj . The branching ratio /i — > 67 is given by 

BR(^-e7) = ^(|A^P + |A^P). (4) 

Owing to the hierarchy of masses in the lepton sector (m^ ^ me) the right-handed term always 
dominates, ^ Hereafter, we will refer to A^^^^ and ^4^"^^ as the chargino and 

neutralino contributions respectively. 
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3 The Universal Case 

In the framework of the constrained MSSM, the occurrence of off diagonal entries in the slepton 
mass matrix m| as introduced in the current setup, gives rise to charged lepton flavor violating 
decays of the type BR{li — > /j7). The corresponding rates depend in general on the supersym- 
metric parameters of the theory, which, we recall, are given on the one side by the GUT-scale 
inputs mo, Mi/2) ^0 and on the other side by tan/3, the sign of (j, and the relative size of the 
relevant off-diagonal entry of which we denote as 5od- In this setting, BR{li — > Ij'^) is, to 

a good accuracy, proportional to tan /?, while the dependence on the sign of /i and on Aq is rather 
weak^. As regards the mass scales of the theory, M1/2 and mo, it is clear that larger masses will 
in general suppress lepton flavor violating decays, since particles circulating into the loops are 
heavier. In particular, being the masses of the charginos and neutralinos mostly insensitive to 
mo, the strongest dependence will be that on Afi/2- 

We plot in the left panel of fig.^the results we get for BR(^ — > 67) as a function of M1/2, in 
the CMSSM scenario for (5od = 10^^, at tan/3 = 10, /i > at various values of tuq. Henceforth, 

^See e.g. the discussion in |12|. 
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Figure 1: : The BR(/i — > 67) {left ) and the chargino-sneutrino {right, sohd red hnes) and 
neutrahno-slepton {right, dashed blue hnes) contributions as functions of the common gaugino 
mass Ml/2 for mo = 250, 500, 750 and 1000 GeV. Aq = 0, tan (3 = 10 and ^ > 0. 

where not otherwise specified, ^0 is set to zero and 5od = 10^'^, for definiteness. We verified that, 
within the present setup, the dependence of LFV branching ratios on Aq is always negligible. 

In the right panel of fig. ^ we split the contributions arising from the chargino and the 
neutralino loops. We observe that, as a rule, the chargino contribution is larger than that of the 
neutralino for sufficiently low values of M1/2, while, vice-versa the neutralino may dominate at 
large Mi/2- Moreover, the chargino contribution is almost independent of mo, on which instead 
the neutralino contribution rather critically depends. This behavior stems from the steepness 
of the different functions entering A^'^^ and A^"'\ whose arguments are the ratios (m-^_/m?) and 

(m^o l^f)-, ill region where they are evaluated. 

Finally, we stress that the two contributions have the same sign, possibly excluding a narrow 
region at very low Mi 12, where the chargino contribution is however always by far dominant. 
Flipping the sign of ^ would change the sign of both contributions simultaneously. Therefore, 
in the universal case no cancellations among the chargino and the neutralino contribution can 
occur. 

4 Non-universal Wino Mass Term M2 

In this section we focus on the effects, on LFV processes, of taking, at the unification scale Mqut, 
the gaugino mass pattern Mi = M3, M2 = r2Mi, and varying r2- The effect of a non- universality 
in the high energy input value of the wino mass term mainly translates, at the low energy scale, 
into a proportional variation of M2, and thus of the wino content of both the charginos and 
the neutralinos ^H]- We remark that, through the electroweak symmetry breaking conditions, 
however, also the value of /u is slightly sensitive to r2. 

We start our analysis studying, in fig. [21 the branching ratio BR(;U 67) as a function of 
r2, at fixed Mi = M3 = 1500 GeV and tan/3 = 10. The sign of is chosen to be positive, 
and the various lines correspond to different choices of mo, ranging from 500 GeV up to 1250 
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Figure 2: The BR(/u 67) as a function of the ratio r2 = M2/M1 for Mi = M3 = 1.5 TeV at 
tan /3 = 10, /X > 0, for mo = 500, 750, 1000 and 1250 GeV. 

GeV. The shaded region around r2 = signals a chargino mass which is not compatible with the 
experimental lower bound th^- > 100 GeV from LEP Jj4j. The extremum right hand-side of the 
plot represents the particular point at r2 = 1 corresponding to gaugino mass universality. We 
notice that no particular effect, a part the obvious suppression of BR(^ 67) due to a heavier 
chargino and neutralino spectrum for increasing \r2\, is present in the region where the sign of 
M2 is the same as the sign of Mi = M3 (i.e. r2 > 0). Contrary, for negative r2, cancellations may 
occur between the chargino and the neutralino contributions. When the two contributions are 
exactly equal and opposite, BR(^ — > 67) drops to zero, as signaled by the sharp dips appearing 
in the figure. In the region between two dips the dominant contribution is that of the neutralino, 
while for r2 ^ and r2 ^ 00 the chargino contribution always dominates. We emphasize that 
such cancellations would not be spoiled by the presence of a subdominant contributing amplitude, 
coming e.g. from Ag, or double mass insertions. The only net effect would be a small shift 
in the precise location of the dips, which would be nonetheless present. 

Fig. 131 shows the separate contributions from the chargino and the neutralino for the two 
extreme cases mo = 500 and 1250 GeV. For the sake of clarity, we also plot the neutralino 
contribution with opposite sign, in order to highlight the points where the cancellations take 
place for mo = 500 GeV. The reason of the wide dip in the mo = 1250 GeV case is also apparent 
from the right panel. 

The location of the dips depends rather sensitively on the value of mo, and in some cases 
it may take place in the cosmologically interesting region where the neutralino is wino-like. In 
some cases (as for mo=1250 GeV in the figure) the exact cancellation never takes place, though 
it is evident that the chargino and neutralino loops interfer destructively, generating a minimum 
at some value of r2, contrary to what one would expect from the overall increase of the mass 
spectrum. 

The chargino contribution depends on the orthogonal matrices Ol and Or which diagonalize 
the chargino mass matrix Mq- In particular, we define Or as the orthogonal matrix which 
diagonalizes McMq, and Ol as that which diagonalizes the combination M^Mc- We call 
chargino masses the two diagonal elements M-- of the matrix OrMqOJ^, which can therefore 



4 



m,=1250GeV 



Figure 3: The separate chargino and neutrahno contributions to BR(// 67) in the two cases 
mo = 500 and 1250 GeV shown in fig. 2. 

also take negative values Multiphcation Ol by the Pauh matrix aUows to recover 

the positive values for the masses. We recall that the diagonal entries of Mc are M2 and 
therefore a large mixing, and sizable off-diagonal entries in Ol^r, appear when M2 = ±/i. Since 
the largest element in the chargino amplitude is given by 

M- 

~ -const • V {Or)^,{Ol)a2^U^iU^2C{xax), = Ml./ml^, (5) 

we expect an enhancement in A^'^^ when M| = being the amplitude always proportional to an 
off diagonal element of Ol,r- Though a peak indeed takes place at M2 = ±/i (around r2 — ±1.4 
in fig. IIJ in the separate contributions from the lighter chargino (A = 1) and from the heavier 
chargino {A = 2), the overall result, shown in the upper panel, is instead completely smooth'^. 
In order to understand why the peaks are smoothed out in the total chargino contribution, a 
closer look at the amplitude © is necessary. The resulting cancellation depends on the relative 
sign of the eigenvalues M-- of Mc and on the products of the matrix elements (0_r)ai(Ol)a2 for 
A = 1,2. We report in tab.^the signs of the various contributions entering eq. ©. Notice that 
the overall contributions have a minus sign in front. For r2 far from zero, the difference between 
the A = 1,2 contributions, which is roughly constant even in the region close to M2 = ±/i, is 
driven by the combination ^^(— l)"^(0/j)^i(OL)^2- While the sign flip is due to the change 
of the sign of M2, the peaks at r2 ~ are due to the sudden decrease of the particle mass 
spectrum. We stress that the absence of a peak at M2 = ±/i is entirely due to the particular 
arrangement of signs shown in tab. ^ An analogous cancellation takes place also in the case of 
the neutralino, though the more involved analytics make it less transparent to elucidate. 

Finally, we comment on the dependence of the outlined structure for lepton flavor violation 
on the other relevant supersymmetric parameters. First, taking a different value for tan/3 
changes the position of the dips, as highlighted in the left panel of fig. El This is due to the 
fact that the chargino contribution is proportional to (1/ cos /3) ~ tan/3, while the neutralino- 
slepton amplitude is proportional to terms of the type (a + b/ cos /3) ^Uj: as tan (3 increases, the 

^Enhancements at M2 = ±/.t have instead been reported in other related frameworks, such as in supersymmetric 
electroweak baryogenesis |16| 
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Table 1: The sign of the terms appearing in the chargino contribution of eq. (4.1). 
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Figure 4: ( Upper panel ): The chargino contribution as a function of r2 at Mi = M3 = 1.5 TeV, 
tan/? = 10, /i > and mo = 1000. {Lower panel ): the contributions from the hght [A = 1) and 
heavy {A = 2) chargino loops at the same supersymmetric parameters as above. Also included 
is the A = 2 contribution with changed sign. 
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neutralino contribution grows therefore more slowly than the chargino contribution. Though 
this effect is rather small, the position of the dips can substantially change, as showed in the 
right part of fig. |S1 

As far as the dependence on the relative size of the off-diagonal entries of m^, is concerned, 
we study in fig.lElthe variations induced in BR(/i ej) by changing, in a range of four orders of 
magnitude, the size of 5oo- We find, as expected, that BR(/x — > 67) oc 6^^^ to a good accuracy. 
Moreover, increasing 6od slightly enhances the chargino contributions with respect to those of the 
neutralino, so that the dips tend to get closer to each other. Nevertheless, no critical dependence 
on 60D is found in the locations of the dips within a span of four orders of magnitude. 

5 Non-universal Gluino Mass Term M3 

Since gluinos do not take part in the loops of LFV decays, the gluino mass term can only 
indirectly affect the process — > 67. Indeed, the electroweak symmetry breaking condition 
gives, at tan (3 = 10, 

+ ^"^1 ~ -0.1ml + 2.1M| - 0.22M| + O.I9M2M3 + O.O3M1M3, (6) 

and the coefficients vary rather mildly over the moderate tan /3 region. We see that the dominant 
contribution to the right-hand side comes from M3; hence, ^ is practically determined by the 
gluino mass. Notice also that for small values of M3, /i^ becomes negative, thus signalling the 
absence of radiative electroweak symmetry breaking. 

In the left panel of fig. [3 we show BR(fj, 67) as a function of r^ for different values of 
mo and /x > 0. In agreement with eq.®, the branching ratio is symmetric with respect to 
r3 ^ — T3. Hence, we can limit our analysis to the case r^ > 0. We notice that the behavior of 
the branching is not monotonous. At small values of r^, BR{fj, — > 67) decreases until it reaches 
a point Tc where it begins to increase. Larger values of mo imply larger values of re- 
in order to understand this peculiar behavior, we have split, in the right panels of fig. El 
the chargino and the neutralino contributions. It is clear from the figure that at small r^ the 
branching ratio is dominated by the chargino contribution, which decreases with r^. At larger 
values of r^ the neutralino contribution dominates, increasing with r^. rc is the point where 
these two contributions become equal. We emphasize that such behavior is rather surprising, 
because at large r^ the particle spectrum becomes heavier, and one would therefore expect, 
instead, some kind of suppression in BR(/x — > 67). 

6 Analysis of BR(yU 57) in Specific Models 

Until now, we have studied generic models with non-universal gaugino masses. In this section, 
we consider two well-motivated specific models which constrain soft-breaking parameters and 
predict explicit relations between the three gaugino masses: SU{5) GUT inspired gaugino non- 
universality and minimal gaugino mediation. 

6.1 SU(5) 

In supergravity models, the gauge kinetic function depends on a chiral superfield (p, whose 
auxiliary F-component acquires a large vacuum expectation value. Gaugino masses come from 
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Figure 5: [Left ): The BR(/i — > 67) as a function of r2 at Mi = M3 = 1.5 TeV, /i > and 
mo = 1000 for tan/? = 3, 10, 30 and 50. [Right ): The details of chargino and neutrahno 
contributions for tan /3 = 3 and 30. The supersymmetric parameters are the same as for the left 
panel. 
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Figure 6: The BR(/x — 67) as a function of r2 at Mi = M3 = 1.5 TeV, > 0, tan/3 = 10 and 
mo = 1000 for various (^od = 10"\ 10"^, IQ-^, IQ-^ 
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Figure 7: (Left ): The BR(|U — > 67) as a function of at various mg. In the figure Mi = 
M2 = 1500 GeV, tan/3 = 10 and ^ > 0. (Right ): The BR(/i ^ 67) as a function of the gluino 
high energy input value M3, parametrized by the quantity = M3/M1 {upper panel ) and the 
separate contributions from the chargino (sohd hue) and neutrahno (dashed hue) loops (lower 
panel ). For definiteness, we set Mi = M2 = 1200 GeV, tan /3 = 10, ^ > and mo = 600 GeV. 



the following five dimensional operator 



L = f> A,A,, (7) 



where Ai^2,3 are the bino B, the wino W and the gluino g fields respectively. In models based 
on SU{5), gauginos belong to the adjoint representation of SU{5), and can belong to any 
irreducible representation appearing in their symmetric product, 

(24 X 24)sym = 1 + 24 + 75 + 200 . (8) 

The mSUGRA model assumes (j) to he a singlet, which implies equal gaugino masses at the GUT 
scale. On the other hand, if cj) belongs to one of the non-singlet representations of S'C/(5), gaugino 
masses are not universal, and they are related to one another via representation invariants [2]. 
The resulting ratios for the gaugino masses at the GUT scale jl7j are shown in Table [2 Some 
phenomenological implications of these models were analyzed, e.g., in ^Sj and [7]. 

We show the behavior of BR(/i ej) as a function of M3 for the different representations 
of SU{5). In this case we do not find any cancellation. For the representations 24 and 200 the 
chargino and neutralino contributions have the same sign, so that they always interfere construc- 
tively. In the representation 75, instead, the signs are opposite, but the chargino dominates over 
the whole parameter space. This, however, does not necessarily imply that such cancellations do 
not take place in a generic SU{5) GUT, because the gauge kinetic function can involve several 
chiral superfields belonging to different representations, and the final gaugino masses would be 
linear combinations of the different contributions. 
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Figure 8: Non universal gaugino masses with high energy ratios lying in the representations of 
the symmetric product of two SU{5) adjoints. In this plot we set rriQ = 1200 GeV, /j, > and 
tan P = 10. 

6.2 Gaugino Mediation 

Models with gaugino mediated symmetry breaking naturally emerge in higher dimensional GUT 
theories ^Hl- In such theories, the MSSM matter fields are localized on a visible brane, while the 
gauge fields can propagate in the bulk of the extra dimensions. Supersymmetry is broken on a 
distant hidden brane, separated from the visible one in the extra dimensions. Extradimensional 
locality forbids direct couplings between the two branes, and therefore suppresses all soft masses 
which involve MSSM matter fields (squark and slepton masses as well as trilinear couplings). 
Gauginos, on the contrary, can directly couple to the source of SUSY breaking, acquiring nonzero 
masses. As a result, below the compactification scale, Mc ~ 1/R, the effective 4-dimensional 
theory is the MSSM, and gaugino masses are the only non-negligible sources of SUSY breaking. 
At low energies, scalar masses and trilinear couplings are generated from gaugino masses through 
renormalization group running. 

Prom the phenomenological point of view, models with gaugino mediation are very attractive. 
On the one hand they alleviate the flavor and CP problems in supersymmetric theories, on the 
other hand, in a context of grand unification, the doublet-triplet splitting of Higgs fields can be 
easily achieved, and proton decay through dimension five operators can be naturally suppressed. 
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Table 2: Relative values of the gaugino masses at the GUT scale for the different representations 
allowed in SU{5) 
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Figure 9: BR(/i 67) in a model with gaugino mediation and non-universal gaugino masses at 
Mc = Mqut- Ml is set to 2 TeV, tan /? = 10 and fi > 0. We plot the results as functions of M2 
at four different values for M3 = 750, 1000, 1250 and 1500 GeV. 
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Figure 10: BR(// ej) in models with gaugino mediation and restricted flipped SU{5) gauge 
symmetry, which enforces non-universal gaugino masses and dictates M2 = M3. 

In these scenarios, at the scale Mc, which we assume to be the unification scale (Mc = Mqut); 
the only free parameters are the gaugino masses Mi and tan /3, since mo and Aq are set to zero. 
In order to avoid a charged LSP (the lightest stau) one has either to resort to some specific GUT 
running above Mc in the universal case [20], or to gaugino non-universality |2HI22[ [23]. In what 
follows, we will limit our discussion to the case, relevant for the present context, of non-universal 
gaugino masses. In figlH we show a plot of BR(// — > 67) as a function of M2 for fixed Mi and 
different values of M3. Once again, cancellations occur when M2 and Afi have opposite signs. 

If SUSY is broken on a brane with restricted gauge symmetry, non universal gaugino masses 
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are generated. For example, if the bulk symmetry is 50(10), hidden branes with Pati-Salam 
SU{4:) X SU{2)l X SU{2)r, Georgi-Glashow SU{5) x U{1), or flipped 5C/(5)' x C/(l)' gauge 
symmetries can be obtained [^H^. If gauginos get masses on these branes, the gauge symmetry 
relates the gaugino masses of the MSSM at the compactification scale. In flipped SU{5), for 
instance, one obtains M3 = M2 7^ Mi [251 123j . In Fig llOl we have considered this case. We 
show the branching as a function of M2 = M3 for several values of Mi. In this specific model, 
for any value of Mi there exists a value of M2 = M3 for which the exact cancellation between 
chargino and neutralino contributions, and the resulting suppression of BR(/i — > ej), occurs. As 
emerging from the figure, this value turns out to be an increasing function of Mi. 

7 Conclusions 

We studied in this paper the implications of non-universal gaugino masses on LFV decays of 
the type li — s- Ij^. The analysis of the interplay between chargino-sneutrino and neutralino- 
slepton loop amplitudes lead us to conclude that cancellations between the two contributions 
may drastically suppress BR(// — > ej) when non-universalities are allowed in the gaugino sector 
at the high energy scale Mqut- 

In particular, we found that such cancellations, which are absent in the case of universality, 
take place when the signs of M2 and Mi are opposite. Moreover, we highlighted that BR(^ 67) 
does not always decrease when the particle mass spectrum increases, but, owing to the neutralino 
contribution dominance, and the above mentioned destructive interferences, it may increase as 
well. In particular, whenMs ^ M2,Mi, BR(/i 67) is found to be an increasmg' function of M3. 
We showed that a peculiar sign arrangement, which we analyzed in some detail, smoothes the 
sharp peak which takes place in the separate contributions from the two charginos at M2 = =b/i. 
We also concentrated on particular models where gaugino non-universality plays a fundamental 
role, as in gaugino mediation, and showed that even in this case where one deals with a restricted 
parameter space, cancellations may greatly suppress BR(;U 67). Further constraints on the 
ratios between the three gaugino masses, stemming from grand unification, may still allow these 
types of suppressions, as we demonstrated in the particular case of flipped SU{5), again in the 
context of gaugino mediation. 

We emphasize that our approach does not depend on the particular model-dependent source 
of lepton flavor violation in the MSSM, as long as the outlined assumptions about the dominant 
source of LFV are fulfilled; the presence of subdominant contributions would not change our 
conclusions. 
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